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Abstract. We present the results of a wide-field survey 
for extremely red objects (EROs hereafter), the widest so 
far, based on Ks and R band imaging. The survey covers 
701 arcmin 2 and it is 85% complete to Ks < 18.8 over the 
whole area and to Ks < 19.2 over 447.5 arcmin 2 . Thanks 
to the wide field covered, a complete sample of about 400 
EROs with R — Ks >5 was selected. The distribution of 
the EROs on the sky is strongly inhomogeneous, being 
characterized by overdensities and large voids. We detect 
at the 8cr level a strong clustering signal of the EROs which 
is about an order of magnitude larger than the clustering 
of if -selected field galaxies in the same magnitude range. 
A smooth trend of increasing clustering amplitude with 
the R — Ks color is observed. These results are strong 
evidence that the largest fraction of EROs is composed 
of high-^ ellipticals, of which we detect for the first time 
the z ^ 1 large scale structure clustering signal. We show 
how the surface density variations of the ERO population 
found in our survey can explain the highly discrepant re- 
sults obtained so far on the density of z > 1 ellipticals, 
and we briefly discuss the main implications of our re- 
sults for the evolution of elliptical galaxies. The number 
counts and the colors of the fGselected field galaxies are 
also presented and briefly discussed. 

Key words: Galaxies: elliptical and lenticular, cD, evo- 
lution, formation, clusters: general, starburst - large-scale 
structure of Universe 
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1. Introduction 

Near-infrared surveys prompted the discovery of a popu- 
lation of objects with very red optical-infrared colors (Ex- 
tremely Red Objects, EROs hereafter; e.g. Elston et al. 
1988, McCarthy et al. 1992, Hu & Ridgway 1994; Thomp- 
son et al. 1999; Yan et al. 2000; Scodeggio & Silva 2000). 
In general, objects have been classified as EROs when they 
had redder colors than late type galaxies (with negligible 
dust extinction) at any redshift. However, depending on 
the depth of the photometry and on the available filters, 
different selection criteria have been used to select EROs. 
In this paper, EROs are defined as objects with R — Ks >5 
(see Sect. 5 for more details on this choice). 

The red colors of EROs are consistent with two classes 
of galaxies: they could be old, passively evolving elliptical 
galaxies at z <L 1 which are so red because of the large 
iv-correction. EROs may also be strongly dust-reddened 
star-forming galaxies or AGN. The observational results 
of the last few years showed that both classes of galax- 
ies are indeed present in the ERO population: on one 
hand, a few objects were spectroscopically confirmed to 
be z > 1 ellipticals (Dunlop et al. 1996, Spinrad et al. 

1997, Liu et al. 2000, and marginally, Soifer et al. 1999), 
or to have surface brightness profiles consistent with be- 
ing dynamically relaxed early type galaxies (e.g. Stiavelli 
et al. 1999, Benitez et al. 1999). On the other hand, other 
EROs have been detected in the sub-mm (Cimatti et al. 

1998, Dey et al. 1999, Smail et al. 1999, Andreani et. al. 
2000), thus providing examples of high-redshift starburst 
galaxies reddened by strong dust extinction and charac- 
terized by high star formation rates. The relative contri- 
bution of the two classes of objects to the whole ERO 
population is still unknown, but there are preliminary in- 
dications, based on near-IR and optical spectroscopy and 
on surface brightness analysis, that ellipticals may repre- 
sent the largest fraction of this population (e.g. Cimatti 
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et al. 1999, 2000; Liu et al. 2000; Moriondo et al. 2000). A 
small fraction of low-mass-stars and brown dwarfs among 
EROs is also expected in case of unresolved objects (e.g. 
Thompson et al. 2000, Cuby et al. 1999). 

The importance of studying EROs is clear especially 
for the clues that they could provide on the formation and 
evolution of elliptical galaxies. For instance, existing real- 
izations of hierarchical models of galaxy formation predict 
a significant decline in the comoving density of the ellipti- 
cals with z, as they should form through merging at z ^ 2 
(Kauffmann 1996, Baugh et al. 1996a), so that a measure 
of a decline of their comoving density would provide a 
stringent proof of these models. Conflicting results have 
been found so far about such issue: some works claim the 
detection of a deficit of z > 1 ellipticals (e.g. Kauffmann 
et al. 1996, Zepf 1997, Franceschini et al. 1998, Barger 
et al. 1999, Menanteau et al. 1999), whereas others hnd 
that a constant comoving density of ellipticals up to z ~ 2 
is consistent with the data (Totani & Yoshii 1997, Ben- 
itez et al. 1999, Broadhurst & Bowens 2000, Schade et 
al. 1999). A potentially serious problem in these studies 
is suspected to be the influence of the field-to-field varia- 
tions in the density of EROs due to the small fields of view 
usually covered in the near-infrared. For instance, Barger 
et al. (1999) found a very low surface density of EROs in 
a 60 arcmin 2 survey to K = 20, while on a similar area 
McCracken et al. (1999) observed a density three times 
larger at the same K level. 

Similar uncertainties have been found in the attempts 
of deriving the fraction of high redshift galaxies among 
IR selected samples, which is believed to be a stringent 
test for the formation of the massive galaxies (Broadhurst 
et al. 1992, Kauffmann & Chariot 1998). Fontana et al. 
(1999), using photometric redshifts, found that the frac- 
tion of high-z galaxies in a collection of small deep fields 
complete to K = 21 was low and comparable with the 
predictions of the cold dark matter (CDM) models, but 
not with passive evolution (PLE) models. The prelimi- 
nary results of Eisenhardt et al. (2000), suggesting a much 
higher fraction of high-z galaxies, are instead consistent 
with both CDM and PLE models. 

The main aim of our survey was to encompass the 
difficulties induced by the cosmic variance, obtaining a 
sample of EROs on a large area, at moderately deep K 
levels, in order to minimize and possibly to detect the 
effects of their clustering, and to compare the observed 
density with that expected in the case of passive evolution 
of ellipticals. Our survey is larger by more than a factor of 
four than the Thompson et al. (1999) survey, and by more 
than an order of magnitude than all the other previous 
surveys for EROs, at the same limiting magnitudes. With 
the large area covered we aimed also to detect a sample, 
or place limits to the surface density of the very rare class 
of extreme EROs with R - Ks > 7. 

In this paper, the observational results of this survey 
and their main implications are presented. A more de- 



tailed interpretation of our findings will be presented in a 
forthcoming paper (Daddi et al., in preparation). 

The paper is organized as follows: we first describe the 
data reduction and analysis, then we present the counts of 
field galaxies. In Sect. [| the sample of EROs is described. 
Sect. H contains the analysis of the clustering of field galax- 
ies and EROs. The main implications of our findings are 
discussed in Sect. ||. H = 50 km s _1 Mpc -1 throughout 
the paper. 

2. Observations, data reduction and photometry 

2.1. Ks-band imaging 

The Ks observations were made with the ESO NTT 3.5m 
telescope in La Silla, during the nights of 27-30 March 
1999, using the SOFI camera (Moorwood et al. 1998) 
with a field of view of about 5'x5'. SOFI is equipped 
with a Hawaii HgCdTe 1024x1024 array, with a scale 
of 0.29"/pixel. The Ks filter has A c = 2.16 [im and 
AA ~ 0.3 \im and it is slightly bluer than the standard K 
filter in order to reduce the thermal background. 

The center of the observed field is at a = 14' l 49 m 29 s 
and 5 = 09°00'00" (J2000). The observed field is one of 
the fields described in Yee et al. (2000) to which we refer 
for details about its selection. The main criteria were not 
to have any apparent nearby clusters and to be at high 
galactic latitude. 

The images were taken with a pattern of fixed offsets 
of 144" (about half of the SOFI field of view) over a grid of 
9x13 pointings. The total area covered by the observations 
was about 24x34 arcmin, with a local integration time 
of 12 minutes in the central deepest region of the field. 
In the shallower region, the effective integration time is 
reduced to not less than 6 minutes. The total amount of 
time required to cover the whole field was about 5.5 hours. 

The data reduction was carried out using the IRAF 
software. The images were flat-fielded with twilight flats. 
The sky background was estimated and subtracted for 
each frame using a clipped average of 6-8 adjacent frames 
(excluding the central frame itself). The photometric cal- 
ibration was achieved each night with the observation of 
5-7 standard stars taken from Persson et al. (1998). The 
zero-points have a scatter of ~ 0.015 magnitudes in each 
night and a night-to-night variation within 0.02 magni- 
tudes. Each frame was scaled to the same photometric 
level correcting for the different zero-points and airmasses. 
Accurate spatial offsets were measured for each frame us- 
ing the area in common with the adjacent frames. The 
images were then combined, masking the known bad pix- 
els, in order to obtain the final mosaic. The cosmic rays 
were detected and replaced by the local median using the 
task cosmicrays of the IRAF package ccdred. The effec- 
tive seeing of the final coadded mosaic ranges from 0.9" 
to 1.1". 
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2.2. R-band imaging 

The i?-band data were taken in May 19-21 1998 with 
the 4.2m William Hcrschel Telescope on La Palma. The 
observations were done using the prime focus camera, 
equipped with a thinned 2048x4096 pixels EEV10 chip, 
with a scale of 0.237"/pixel. This gives a field of view of 
about 8.1' x 16.2'. A standard Johnson i?-filter was used. 
The whole field has been covered by a mosaic of 6 point- 
ings. Each pointing consisted of at least 3 exposures of 
1200s taken with small offsets. The total integration times 
per pointing was therefore 3600s, with the exceptions of 
two pointings with 4800s and 6000s. 

The photometric calibration was achieved with stan- 
dard stars taken from Landolt (1992) with a scatter in the 
zeropoints, from the different stars used, below 0.01 mag- 
nitudes. The images were de-biased and then flatfielded, 
using a master flatfield constructed from the science ex- 
posures, scaled to the same zeropoint and then combined. 
The seeing of the final i?-band mosaic was between 0.7" 
and 0.8". 

2.3. Sample selection and Ks photometry 

The software SExtractor (Bcrtin & Arnouts 1996) was run 
on the Ks mosaic with a background weighted threshold 
in order to take into account the depth variations across 
the defined by SExtractor. Among the detected 

objects, all those with S/N> 5 in a 2" circular aperture 
(twice the average seeing FWHM) were selected and added 
to the catalog. A few spurious detections (e.g. close to im- 
age defects) have been excluded after a visual inspection 
of the image. The final catalog includes 4585 objects. In 
the central deepest region, the 5a limiting aperture magni- 
tude is _Ks(2")=19.6, whereas in the remaining area such 
limit is Ks(2") £ 19.2 because of the reduced integration 
time. 

Isophotal magnitudes were measured with a limiting 
threshold of about 0.7cr s k y corresponding in the central 
area to a surface brightness limit of about ~ 21 mag 
arcsec™ 2 . The aperture correction from 2" to total mag- 
nitudes was estimated throughout the area by measuring 
the difference between the isophotal and the 2" aperture 
magnitudes for the stars with Ks <16. A differential cor- 
rection, in the range of 0.16-0.30 magnitudes, was mea- 
sured for different regions of the mosaic with a typical 
scatter less than 0.03 magnitudes. For the bright objects 
the isophotal magnitudes were on average consistent with 
the Kron automatic aperture magnitudes. However, we 
adopted the isophotal magnitudes because the Kron mag- 
nitudes are rather unstable at faint flux levels, where the 
low signal often does not allow to define the correct auto- 
matic aperture. 

The total Ks magnitudes were then defined as the 
brightest between the isophotal and the corrected aperture 
magnitude. This allowed to safely assign a total magni- 



tude for both the faint and the bright objects. The typical 
Ks magnitude where the corrected aperture magnitude 
begins to be adopted as the total one is Ks~18 in the 
central deepest region. 

The completeness of our catalog has been estimated by 
adding artificial objects to the Ks mosaic in empty posi- 
tions, using the IRAF package artdata. Point-like sources 
as well as objects with De Vaucouleurs and exponential 
profiles (convolved with the seeing PSF) were simulated, 
and SExtractor was run with the same detection parame- 
ters as for the real data. The 85% completeness magnitude 
for the deepest area is Ks = 19.2 for point-like sources. 
The completeness decreases to ~ 70% for the worst case 
that we have tested, i.e. for exponential galaxies with 0.7" 
half light radius. In the shallower area the corresponding 
limiting magnitude is Ks < 18.8. Most of the Ks ^ 18.5 
galaxies are anyway expected to be only barely resolved 
with the Ks 1" seeing (Saracco et al. 1997), and this cer- 
tainly occurs for the distant z £ 1 ellipticals, and thus 
their completeness limits can be assumed to be similar to 
those for stars. 

2.4- R-band photometry and colors 

In order to recover the i?-band counterparts of the Ks- 
selected objects, a coordinate mapping between the Ks 
and the R images was derived. SExtractor was then run 
in ASSOC mode with a search box of 2 x FWHM R . The re- 
gions around bright stars or defects in the R and Ks band 
images were excluded from this analysis. The final effec- 
tive area is 701 arcmin 2 at Ks < 18.8 and 447.5 arcmin 2 
at 18.8 < K s < 19.2. 

Whenever an object had S/N<3 in the i?-band image 
3cr limits were assigned. The 3cr limiting magnitude in a 
2" diameter aperture is R > 26.2 for most of the area, 
reaching R > 26.5 in the deepest pointing. When S/N>3, 
2" diameter corrected aperture magnitudes were assigned 
to each object. The aperture correction was derived in the 
same way as for the ifs-band, with slightly smaller correc- 
tions because of the better i?-band seeing. The magnitudes 
were dereddened for Galactic extinction. At the Galactic 
coordinates of the center of our field (/ ~ 5?5, b ~ 57°), the 
extinction coefficients from Burstein & Heiles (1982) and 
from Schlegel et al. (1998) are A B = 0.04 and A B = 0.13 
respectively. Since the two values are derived in differ- 
ent ways, neither of the two can be discarded. The aver- 
age was therefore adopted, obtaining a correction of 0.052 
magnitudes in R and negligible in Ks. This introduces 
an uncertainty of ~ 0.03 magnitudes in the dereddened R 
magnitudes. 

Finally, the R — Ks colors are defined for all the ob- 
jects as the difference between the R and Ks corrected 
aperture magnitudes. Thanks to the depth of the .R-band 
data, colors as red as R — Ks = 7 could be measured 
down to the Ks magnitude limits of our survey. Because 
of the long integrations used in the i?-band the objects 
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Fig. 1. The observed number counts for both stars and galax- 
ies in our survey. The error bars indicate the poissonian uncer- 
tainties. 



with R £ 20 were saturated. This effect is obviously more 
important for stars than for galaxies; moreover, since we 
are interested in the extremely red galaxy population, the 
effects due to saturation have no impact on our results. 

2.5. Star-Galaxy classification 

The star-galaxy (S/G) separation was done by means of 
the SExtractor CLASS_STAR parameter in both the R 
and Ks band. This classification was found to be reli- 
able for objects with Ks ^ 17.5 and R ^ 24. Because of 
the seeing variations through the area in both bands, a 
variable CLASS_STAR threshold was adopted in different 
subareas. Given the better seeing of the i?-band, the clas- 
sification was based mostly on that band, switching to the 
Ks CLASS_STAR for objects close to the saturation level 
in R. Objects which are fainter than Ks = 17.5 and R = 
24 have not been classified and have been considered to be 
galaxies. This means that a S/G separation can be pro- 
vided only for objects with colors R — Ks ^ 5 at K s=19, 
R — Ks H 6 at Ks=l8, and so on (see Fig. || and the upper 
panel of Fig. ||). From the small number of objects classi- 
fied as stars which are close to the diagonal straight line 
which indicates our S/G classification limit in the color - 
magnitude plane, we can safely conclude that our inability 
to properly classify very red, faint objects has almost no 
effect on the total star number counts and, therefore, on 
the galaxy number counts. In addition, it is possible that 
near the faint limit of our survey a very small fraction of 
very compact, blue galaxies, such as for instance AGN or 



compact narrow emission line galaxies (e.g. Koo & Kron 
1988), could have been incorrectly classified as stars. 

3. K s-band number counts 

Galaxy number counts in the Ks band can provide more 
advantages in studying galaxy evolution and cosmological 
geometry than optical counts because they are much less 
sensitive to the evolution of stellar population and to the 
dust extinction. Our survey, which covers the magnitude 
range 14 < Ks < 19.2, represents the widest among the 
previous deep surveys at levels fainter than Ks > 18. 

Table 1. Differential Number Counts 



Ks 


range 


Area' 2 ' 


Galaxies 


Stars 


11.7 


- 12.2 


701 




4 


12.2 


- 12.7 


Yl 




7 


12.7 


- 13.2 


y> 




9 


13.2 


- 13.7 


y> 




17 


13.7 


- 14.2 


y> 


4 


21 


14.2 


- 14.7 


Yl 


6 


38 


14.7 


- 15.2 


Yl 


16 


37 


15.2 


- 15.7 


yi 


30 


62 


15.7 


- 16.2 


yi 


74 


73 


16.2 


- 16.7 


yi 


100 


88 


16.7 


- 17.2 


Yl 


178 


128 


17.2 


- 17.7 


Yl 


372 


127 


17.7 


- 18.2 


Yl 


633 


144 


18.2 


- 18.7 


Yi 


892 


185 


18.7 


- 18.8 


Y) 


200 


32 


18.8 - 


- 19.2W 


447.5 


628 


84 



The last bin includes only the objects in the deepest region. 



Table summarizes the number of galaxies and stars 
detected in each Ks bin and Fig. [i] shows the correspond- 
ing differential counts in 0.5 magnitude bins. No correction 
for incompleteness was applied. Galaxies start to domi- 
nate over stars at Ks ~ 16.5 and their surface density is 
about a factor of 8 higher than the stellar surface density 
at Ks ~ 19. 

The slopes of the galaxy number counts were derived 
over the magnitude range covered by our survey. At bright 
magnitudes a slope of 7 ~ 0.53 ±0.02 is found in the range 
14 < Ks < 17.5. We confirm that the if-band galaxy 
counts show a flattening at Ks ~ 17.5, where the best 
fit slope changes from 7 = 0.53 to 7 = 0.32 ± 0.02 (see 
Fig. [2]). The leveling off of the counts below a slope of 
0.4 indicates that the differential contribution to the ex- 
tragalactic background light (EBL) in the K band peaks 
at Ks ~ 17—18 and then starts to decrease at fainter 
fluxes. The contribution to the EBL over the magnitude 
range 14 < Ks < 19.2 sampled by our survey is about 4.20 
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Fig. 2. Counts of galaxies from our survey compared to a col- 
lection of published data. The solid curve shows the average 
counts as estimated by Hall & Green (1998). 

nW/m 2 /sr, which constitutes about 53% of the estimated 
EBL from discrete sources in the if-band (cf. Pozzetti et 
al. 1998, Madau & Pozzetti 2000). 

Fig. |] shows the differential galaxy number counts in 
our survey compared with a compilation of X-band pub- 
lished surveys. No attempt was made to correct for dif- 
ferent filters. As shown in the figure, our counts are in 
very good agreement with the average counts of previous 
surveys (Hall & Green 1998). 



Table 2. Galaxy Median Colors 



Ks 


range 


Galaxies 


Median R - Ks 


16.5 


- 17.0 


143 


3.54 


17.0 


- 17.5 


277 


3.80 


17.5 


- 18.0 


522 


3.92 


18.0 


- 18.5 


759 


4.08 


18.5 


- 18.8 


613 


4.11 


18.8 


- 19.2 


628 


4.04 



4. The sample of EROs 

In Fig. H the R — Ks vs. Ks color - magnitude diagram is 
plotted for both stars and galaxies in our sample. The di- 
agonal straight line in the upper panel indicates our S/G 
classification limit in this plane. Because objects above 
this line are not classified, no star can appear in the up- 



Fig. 3. Color magnitude diagram for the stars (upper panel) 
and galaxies (lower panel) samples. The diagonal straight line 
in the upper panel indicates our S/G classification limit in the 
color - magnitude plane (see text). The bluer objects plotted 
with filled squares have at least one pixel close to saturation in 
the 7?-band. Some of the bluest saturated galaxies at Ks > 17 
may actually be stars. The reddest galaxies with R — Ks ^ 7 
have color lower limits as they are not detected in the i?-band 
(see Fig. . 



per right corner. However, this figure shows that there are 
very few stars with a color redder than R — Ks > 5, even 
in the region of this color - magnitude diagram where our 
morphological classification is still reliable. This suggests 
that very few stars should be present in the sample of 
objects for which no morphological classification was pos- 
sible. 

Fig. U shows the distribution of colors for galaxies 
fainter than Ks — 16.5, from which the median color at 
different Ks magnitude levels have been calculated (see 
Table ||). In spite of the presence of objects with color 
upper and lower limits, the use of the median colors (in- 
stead of the mean colors) allows unbiased estimates in 
the Ks range we are considering. The faintest bin with 
18. 8< Ks <19.2 includes only the objects detected in the 
deeper region. The median R — Ks color of the galaxies 
increases by 0.5 magnitudes from Ks = 16.5 to Ks ~ 18 
and then it remains almost constant up to the limits of our 
survey (Ks=19.2). This trend is similar to what is found 
by Saraceo et al. (1999) for the median J — Ks galax- 
ies color, which also reaches a maximum at Ks ~ 18 19 
and then it becomes bluer, while the median B — K color 
gets significantly bluer at brighter (K ~ 17) magnitudes 
(Gardner et al. 1993). 
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Table 3. The Sample of EROs 



area R - Ks>5 R - ,Ks>5.3 R - Ks>6 R - Ks>7 

Ks limits arcmin 2 N Frac. Dens. N Frac. Dens. N Frac. Dens. N Frac. Dens. 



Ks <17.0 


701 


2 


0.006 


0.003 




















Ks <17.5 


701 


15 


0.025 


0.02 


5 


0.008 


0.007 














Ks <18.0 


701 


58 


0.051 


0.08 


19 


0.017 


0.027 














Ks <18.5 


701 


158 


0.084 


0.23 


75 


0.040 


0.11 


7 


0.004 


0.01 








Ks <18.8 


701 


279 


0.111 


0.40 


150 


0.060 


0.21 


30 


0.012 


0.04 


2 


0.0008 


0.003 


Ks <19.0 


447.5 


220 


0.116 


0.49 


133 


0.070 


0.30 


33 


0.017 


0.07 


4 


0.0021 


0.009 


Ks <19.2 


447.5 


281 


0.126 


0.63 


173 


0.079 


0.39 


44 


0.020 


0.10 


5 


0.0023 


0.011 


Ks <19.2 


701+447.5 


279+119 


0.127 


0.67 


150+87 


0.076 


0.40 


30+27 


0.018 


0.10 


2+4 


0.0019 


0.012 



We present in detail the cumulative number (N) of EROs selected at each Ks limiting magnitude, the fraction of EROs with 
respect to the whole held galaxies (Frac.) and the corresponding surface density (Dens., in objects/arcmin 2 ). The last line was 
calculated using the whole survey area to Ks = 18.8 and the deeper area to Ks — 19.2. The data presented here are used 
throughout the paper. 



17.0 - 17.5 



17.5 - 18.0 



nJW Ik srff 



18.5 - 18.8 



18.8 - 19.2 



0246024 6 
R-Ks 

Fig. 4. Distributions of the R — Ks colors for the galaxies of 
our sample. 



Our wide-field survey allows us to select a statistically 
significant and complete sample of EROs which can pro- 
vide stringent constraints on the density of high-z ellip- 
ticals (see Sect. 1). For this reason, EROs are defined as 
objects with R — Ks > 5 because this corresponds to se- 
lect passively evolving ellipticals at z £ 0.9 (see Fig. |i"3| ). 
Table || shows the results of the selection for different 
magnitude and color thresholds. Among the others, the 
threshold R — Ks > 5.3 is used because it corresponds 
to the selection of z ^ 1 elliptical galaxies (see Fig. |i~3|) . 
We detected 279 EROs with Ks < 18.8 from the whole 
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Fig. 5. Enlarged portion of Fig. |^ (bottom) around the reddest 
colors. Filled boxes are objects within the completeness limits 
of our survey (Ks < 19.2 in the deeper area and Ks <18.8 out- 
side) while filled triangles are objects fainter than those limits. 
The horizontal long-dashed lines correspond to the limits re- 
ported in Table |^ for the selection of the samples of EROs. In 
the region above the short-dashed line the star-galaxies sep- 
aration is not feasible. The objects along the diagonal dotted 
lines are not detected in R and have a 3<j limit in that band. 



area and 119 EROs with 18.8 < K s < 19.2 in the deeper 
area, yielding a total sample of 398 objects (see Fig. ||and 
Table |3|). This is by far the largest sample of EROs ob- 
tained to date. A small but complete sample of EROs with 
R — Ks > 7 has also been selected, and we estimate for 
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the first time their surface density to be ~ 0.01 arcmin 2 
at Ks ~ 19. 

A comparison of the surface densities of the EROs in 
our sample with those in the Thompson et al. (1999) sur- 
vey can be directly done after taking into account the 
different filters used in the two surveys. For the K' filter 
used by Thompson et al. (1999), Ks ~ K' - 0.2 (adopting 
H — K = 1), and therefore their limit at K' < 19 corre- 
sponds to Ks < 18.8 which is the shallower limit of our 
survey, while their R — K' color is bluer than our R — Ks 
by about 0.1 magnitudes for the redder objects (Thomp- 
son, private communication). At the level of Ks < 18.8 we 
find a density of 0.042 ± 0.008 arcmin" 2 for EROs with 
R — Ks > 6, to be compared with the value of 0.039±0.016 
that they find (these errors are poissonian). Thus, the 
two surface densities are in excellent agreement with each 
other. We also verified that the average R — Ks color of all 
our objects with 17.8 < K s < 18.8 (R- Ks = 3.70 ± 0.03, 
determined with a Kaplan-Meier estimator), is in good 
agreement with the average R — K' = 3.73 ± 0.04 in 
18 < K' < 19 by Thompson et al. (2000). 

5. The angular correlation functions 

Statistical measurements of the clustering of faint galax- 
ies are important for studying the evolution of galaxies 
and the formation of structures in the Universe. In fact, 
the amplitude of clustering in 2D space is a useful probe 
of the underlying 3D structure (e.g. Connolly et al. 1998, 
Efstathiou et al. 1991, Magliocchetti & Maddox 1999). 
The clustering of galaxies on the sky has been studied 
extensively especially in the optical, but also in the near- 
infrared (e.g. Roche et al. 1998 and 1999, Postman et al. 
1998, Baugh et al. 1996b). Our survey, as noted before, is 
the widest at the limits of Ks ~ 19. It is therefore inter- 
esting to estimate the clustering of our sample of galaxies. 

5.1. Calculation technique 

The angular two-point correlation function w{9) is defined 
as the excess probability (over a poissonian distribution) 
of finding galaxies separated by the apparent distance 9: 



dp = n 2 [\ + w(e)]dn 1 dn 2 



(i) 



where N is the mean density per steradian (Groth & Pee- 
bles, 1977). 

Several methods for estimating w(0) from a set of ob- 
ject positions have been proposed and used, but the most 
bias-free and suitable for faint galaxies samples resulted to 
be the Landy & Szalay technique (Landy & Szalay 1993, 
see also Kerscher et al. 2000). This technique (adopted 
for the calculations in this paper) consists in deriving the 
counts of objects binned in logarithmic distance intervals, 
for the data-data sample [DD], the data-random sample 
[DR] and the random-random sample [RR] . These counts 
have to be normalized, i.e. divided for the total number 



of couples in each of the 3 samples. From them we can 
estimate w b {9) as: 



w h (d) 



[DD] - 2[DR] + [RR] 
[RR] 



(2) 



which is biased to lower values with respect to the real 
correlation function w(9): 



w{9) = w h (9)+a 2 



(3) 



where a 2 is the "integral constraint" (Groth & Peebles, 
1977): 



<r 2 = i / / inffhiihdih 



(4) 



Assuming that the angular correlation function w{9) can 
be described by a power law of the form w{9) = A9~ s , 
then, following Roche et al. (1999), we can extimate the 
ratio between a 2 and the amplitude A using the random- 
random sample: 



a 2 ENri^Ll 

A E^rr(0) 



(5) 



The amplitude of the real two-point correlation func- 
tion w(9) can then be estimated by fitting to the measured 
Wb(6) the function: 

w b (9) = A(9- 5 



C) 



(6) 



The errors can be estimated, following Baugh et al. 
(1996b), as: 



Sw h (9) = 2 v / (l + w b (9))/DD) 



(7) 



where DD is the non normalized histogram of [DD]. Eq. 
(Q) is equivalent to assuming 2a poissonian errors for the 
correlations, and it gives estimates that are comparable to 
the errors obtained with the bootstrap technique (Ling, 
Frenk & Barrow 1986). This is necessary because it is 
known that, as the counts in the different bins are not 
completely independent, assuming the la poissonian er- 
rors would result in an underestimate of the true variance 
of the global parameters of the angular correlation (see 
Mo et al. 1992). 

In case of the presence of a randomly distributed spuri- 
ous component among the analyzed sample of objects (an 
example of this case would be a residual stellar compo- 
nent among the galaxy sample), the resulting amplitudes 
are apparently reduced by a factor (1 — /) 2 , where / is 
the fraction of the randomly distributed component (see 
e.g. Roche et al. 1999), and the corresponding correction 
should be applied. 

The random samples used in our analysis were ob- 
tained using the pseudo-random number generator rou- 
tine of the C Library function drand^8. Random samples 
with up to 200 000 objects were used. Typically the num- 
ber of objects in the random samples were a factor of 
100-200 larger than the number of observed objects. The 



8 



Daddi et al.: Clustering of EROs 



- 

_ 










Ks S19.2 
» ^ • 


0.1 
0.01 


- ■ 




• 














1 






i — i i i 1 1 1 1 

• 


— i — i — 


1 INI 1 — Hr- 

Ks ^18.8 


0.1 






• 






0.01 








t 




1 

0.1 






1 1 1 1 1 1 1 1 

• 


— i — i — 


Nil 1 — H- 

Ks S18.5 


0.01 












1 

0.1 






i — i i i 1 1 1 1 
• 


— i — i — 


Mil 1 — H- 

Ks 518.0 


0.01 







0.001 



0.01 
0(deg) 



0.1 



Fig. 6. The observed, bias-corrected two-point correlation 
functions of our if-selected sample. The lines plotted are the 
best-fitted power laws (see Table aj. The error in each bin is 
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two times the poissonian error (see Sect. 5.1) 



Fig. 7. The clustering amplitudes measured in our survey, 
compared with published data from the literature. The models 
shown here are from Roche et al. (1998) and Fig. 7 of Roche 
et al. (1999). 



Table 4. Clustering amplitudes for the /^-selected sample 



Ks limit 


area 


Galaxies 


A[HT 3 ] 


S 


C 


18.0 


701 


1131 


1.3±0.5 


0.8 


4.55 


18.5 


701 


1890 


1.6±0.3 


0.8 


4.55 


18.8 


701 


2503 


1.5±0.2 


0.8 


4.55 


19.2 


447.5 


2222 


1.6±0.2 


0.8 


5.16 



random sample was generated with the same geometrical 
constraints as the data sample, avoiding for instance to 
place objects in the regions excluded around the brightest 
stars. 

5.2. The clustering of the K -selected field galaxies 

In our analysis a fixed slope of S = 0.8 was assumed, as 
this is consistent with the typical slopes measured in both 
faint and bright surveys (e.g. Baugh et al. 1996b, Roche 
et al. 1996, Maddox et al. 1990), and because it gives us 
the possibility to directly compare our results with the 
published ones that are typically obtained adopting such 
a slope. The factor C was estimated (with Eq. for 
both the whole and the deeper areas, turning out to be 
4.55 and 5.16 respectively (the angles are expressed in de- 
grees, if not differently stated). In Fig. [|the observed, bias 
corrected, two-point correlation functions w(6) are shown; 



the bins have a constant logarithmic width (A log 9 = 0.2) , 
with the bin centers ranging from 3.6" to 15'. 

We clearly detect a positive correlation signal for our 
sample with an angular dependence broadly consistent 
with the adopted slope S = 0.8, even if the measurements 
show some deviations, in particular for the brightest sam- 
ples. A few cluster candidates are present in our survey. 
These possible clusters include galaxies with R—Ks < 4.5, 
and are therefore expected to be at z ^ 0.6. A detailed 
analysis of the cluster candidates will be given in a forth- 
coming paper. For the purpose of the present work, we 
tested that the measured clustering amplitudes are stable 
in case of removal of the galaxies of the most evident clus- 
ter from the sample. However, the presence of such clus- 
ters, most of which happen to be in the shallower area, 
could partly explain the observed deviations from the fit- 
ted w(9) = A 6~ 8 power laws for the three brightest 
samples. 

The derived clustering amplitudes are presented in Ta- 
ble ||. The amplitude errors are obtained from the fit as- 
suming Eq. (|7|). No correction for the stellar contamina- 
tion was applied. In Fig. |?] the clustering amplitudes of 
our samples are compared with other published measure- 
ments. The data are shown together with a number of 
PLE models with different clustering evolutions, which are 
described in detail in Roche et al. (1998, 1999). Our mea- 
surements are in good agreement with both the models 
and the previous estimates of Roche et al. (1999), except 
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Fig. 8. The sky positions of the objects with Ks <18.8 and 
R — Ks>5 in our survey. The region enclosed within the dashed 
line is the deeper region. The colors of the objects inside the 
circle are shown in Fig. 



for the point with limiting magnitude Ks = 18.0, which 
is however the most uncertain of our data points. 

As a check it was verified that the correlations of the 
stellar sample are consistent with zero, within the mea- 
sured errors, at all the scales. This is a confirmation that 
the stars are homogeneously distributed on the field (as 
they should be) and the seeing variations across the area 
did not cause a detectable bias in our classification. 

5.3. The clustering of the extremely red objects 

The large sample of EROs derived from our survey allowed 
us for the first time to estimate their clustering properties. 
Even a simple visual inspection of the sky distribution of 
the objects with R — Ks > 5 (see Fig. ||) shows that the 
EROs have a very inhomogeneous distribution. 

The results of the quantitative analysis of the cluster- 
ing are shown in Fig. ^ where the observed, bias-corrected 
angular correlations w{9) of the objects with R — Ks>5 
are displayed. A strong clustering is indeed present at all 
the scales that could be measured, and its amplitudes (Ta- 
ble ^|) are about an order of magnitude higher than the 
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Fig. 9. The observed, bias-corrected two-point correlations 
for the sample of EROs (with R — Ks>5) in our survey. As 
in Fig. ^, the error in each bin is twice the poissonian one. Be- 
cause of the small number of objects included, some bins for 
the two brightest samples, at small angular separation, were 
not populated. The corresponding upper limits would be in 
agreement with the fitted amplitudes. 



ones of the field population at the same Ks limits. The 
correlations are very well fitted by a S = 0.8 power law. 
No attempt was made to correc t th e amplitudes for the 
stellar contamination (see Sect. 5.1), and we stress that 
such corrections would increase them. Adopting the er- 
rors derived from the fits, our detections are significant at 
more than 7a level for the samples with R — Ks > 5 and 
Ks limits equal to or fainter than Ks = 18.5. 

The amplitudes shown in Fig. 5 suggest a possible 
trend of decreasing strength of the clustering for fainter 
EROs: the Ks <19.2 EROs are less clustered than the 
ones with Ks <18.5 and the difference is significant at 
2.7cr, based on the derived errors. The significance of this 
effect would, however, decrease if the contamination by 
randomly distributed stars increases towards the limit of 
our survey. Although difficult to be quantified, this is likely 
to happen because, differently from the brightest EROs, 
only a small fraction of the EROs fainter than Ks = 18.5 
could be morphologically classified (see Fig. ||). 

Defining redder thresholds drastically reduces the 
number of EROs and it is not possible to estimate with 
sufficient accuracy how the amplitudes change for objects 
with even redder R — Ks colors. We could only verify that 
the sample of EROs with R — Ks > 5.3 has clustering am- 
plitudes consistent with those of the R — Ks > 5 samples 
(see Table |). To measure the clustering of the R — Ks >6 
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Table 5. Clustering amplitudes for the Extremely Red Objects 





area 


R — Ks > 5 sample 


R — Ks > 5.3 sample 






Ks limit 


arcmin 2 


Galaxies 


A[10" 3 ] 


Galaxies 


A[10~ 3 ] 


5 


C 


18.0 


701 


58 


24±10 






0.8 


4.55 


18.25 


701 


106 


25±5 






0.8 


4.55 


18.5 


701 


158 


22±3 






0.8 


4.55 


18.8 


701 


279 


14±2 


150 


14±3.4 


0.8 


4.55 


19.2 


447.5 


281 


13±1.5 


173 


12±2.3 


0.8 


5.16 



EROs, an area at least 10 times larger than ours (i.e. ^2 
square degrees) at K = 19 would be needed, assuming 
that their clustering amplitudes are similar to those of 
the EROs with R - K s > 5. 

Finally, it was studied if and how the clustering am- 
plitude changes as a function of R — Ks for the Ks <18.8 
sample (see Fig. |l0]). A clear increase of A with R — Ks is 
present for colors R— K s>3.5, while the R—Ks >3 sample 
has an amplitude that is consistent with that of the whole 
sample of field galaxies. The variation of A can be de- 
scribed with a power law in the range of 3 < R—Ks < 5.7. 
Previous efforts to disentangle the clustering properties of 
the red and blue populations in faint iC-selected samples 
probably failed because the ERO population was not suf- 
ficiently sampled. For instance, Kummel & Wagner (2000) 
did not find significant differences in the clustering of ob- 
jects with color bluer or redder than R — Ks=3A9 for their 
K <17 sample. This is not surprising since at K<17 the 
ERO population is almost absent (see Table [| and Fig. |J). 

To check for the stability of these results, possible sys- 
tematics that could produce a bias in our work were an- 
alyzed. First af all, as the clustering of our i^s-selected 
galaxies is in good agreement with the literature data (Fig. 
0) , we can exclude the presence of measurable biases com- 
ing from the selection of the sample. 

Regarding the color measuraments, since EROs are the 
tail of objects in the R — Ks color distribution, systematic 
variations of the photometric zeropoints across the area 
could have the effect of creating artificial ERO overdensi- 
ties and voids. To exclude this possibility we verified that 
the blue tail of the R — Ks distribution is homogeneously 
spread across our survey, with a very low clustering ampli- 
tude. In case of zeropoints variations these should produce 
the same effect in both the tails of the color distribution. 
Moreover, to test the reality of the large void of EROs 
clearly seen in the bottom region of our survey (see Fig. 
||), the R — Ks color distribution of the galaxies inside 
and outside this large void were compared by means of a 
Kolmogorov-Smirnov test, selecting only the galaxies with 
R — Ks <4 in both regions. The probability that the two 
distributions are extracted from the same population is 
43%. Thus, the two regions are fully consistent with each 
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Fig. 10. The measured correlation amplitude for the samples 
with R — Ks redder than a fixed threshold, in function of the 
threshold itself. This was obtained with the Ks <18.8 sample. 
The horizontal lines show the ±lcr range for the clustering of 
the whole population of field galaxies at this level. The fitted 
power law (dashed line) is log A = 0.43(7? - Ks) - 4.04. 



other with respect to the color distribution of the blue 
population. This, together with the fact that no underden- 
sity is present in this region when only the bluer galaxies 
are considered, shows that the void of EROs should be 
considered a real feature. 

All these tests strongly suggest that the inhomoge- 
neous ERO distribution is a real effect. 



6. Main implications 

6.1. On the nature of EROs 

The strong clustering signal that we find to increase with 
the R — Ks threshold and to reach very high values for 
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the EROs is potentially capable of giving insight about 
the nature of these objects. 

The main possible source of objects that may con- 
tribute to the ERO population, as discussed in the intro- 
duction, are old passively evolving z £ 1 ellipticals, dust- 
reddened starburst galaxies and, in case of unresolved ob- 
jects, low-mass stars or brown dwarfs. Our field being at 
high galactic latitude (b ~ 57°), stars are expected to have 
no clustering, and to be homogeneously distributed, and 
certainly not to give the strong clustering signal detected. 
As for the starburst galaxies, it must be noted that in such 
galaxies the red colors are mainly driven by the amount 
of dust extinction and not by the redshift, as in the case 
of ellipticals (see Fig. [l3]), and therefore a wide redshift 
distribution is expected which should dilute their intrinsic 
clustering. Moreover, it is known that the IRAS-selected 
galaxies (which are typically star-forming galaxies) have 
very low intrinsic clustering (e.g. Fisher et al. 1994). We 
can therefore reasonably conclude that the observed signal 
is due to the clustering of high redshift ellipticals. This is 
also suggested by studies of the local universe which have 
shown that early-type galaxies are much more clustered 
than late-type galaxies (e.g. Guzzo et al. 1997, Willmer et 
al. 1998). In this regard the results plotted in Fig. |l0| could 
be qualitatively explained by noting that in selecting red- 
der samples the fraction of early type galaxies increases 
(the color of local ellipticals is just around R — Ks ;$ 3) 
and by assuming that these galaxies are intrinsically more 
clustered, while such plot would be difficult to understand 
if mainly driven by the strongly reddened starburst galax- 
ies. These considerations strongly suggest that EROs are 
mainly composed by z ^ 1 ellipticals, confirming the pre- 
vious indications that had been found on this issue. 

As the elliptical galaxies are the dominant population 
of galaxy clusters, we investigated the possibility that the 
detected clustering of EROs could be the result of a few 
massive clusters at z > 1 present in our field. For ex- 
ample, in the region inside the circle in Fig. [§, a large 
ERO overdensity is found, that one could suspect to be 
due to a high-z cluster of galaxies. However, Fig. [Tl| shows 
that there is no clear color-magnitude sequence among the 
R—Ks > 5 objects inside that region, suggesting that they 
do not all belong to a single cluster. In case of a cluster, 
even at high-z, a well defined color-magnitude sequence 
is in fact generally observed (e.g. Stanford et al. 1998). In 
the last years a few examples of massive z ^ 1 clusters of 
galaxies have been discovered (e.g. Stanford et al. 1997, 
Rosati et al. 1999), with X ray luminosity of ~ 10 44 erg 
s _1 . A crude estimate of the number of structures of this 
sort that could be observed in our survey can be derived by 
calculating the number of high-z clusters with Lx > 10 44 
erg s _1 expected in the volume we are sampling. From 
the X ray luminosity function of such structures at z ~ 1 
(Rosati et al. 2000) we estimate that the expected num- 
ber of massive clusters in our field in the redshift interval 
0.9 < z < 2 is only ~ 0.1 (for fl = 1). 
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Fig. 11. Color magnitude diagram of all the objects with 
Ks <18.8 in the region inside the circle of Fig. ^. 

Moreover, the detection of the ERO positive correla- 
tion, following a 8 — 0.8 power law on all the scales from 
10" to 15' (corresponding to ^8 Mpc at Z ~ 1) suggests 
that the clustering signal does not come from a few possi- 
ble clusters detected in our field, but rather from the whole 
large scale structure traced by the elliptical galaxies. 

6.2. Fluctuations of the ERO number density 

Our results on the clustering of EROs have important 
consequences on the problem of estimating the density 
of high-z ellipticals (see Sect. 1). 

The existence of an ERO angular correlation with 
8 = 0.8 and with a high amplitude implies significant sur- 
face density variations around the mean value even for rel- 
atively large areas. In the presence of a correlation with 
amplitude A, the rms fluctuations of the counts around 
the mean value n is (see for example Roche et al. 1999): 



n (1 + nAC) 



(8) 



The factor C is the same as in Eq. (|J) and, by apply- 
ing Eq. (||) for several areas, it was found that it can be 
approximated as: 



C = 58 Area 



-0.4 



(9) 



if the area is expressed in arcmin 2 and 8 = 0.8. The valid- 
ity of such an approximation has been tested for square 
regions and for areas not larger than the ones of our sur- 
vey. With Eq. (||) and (^) the expected variations of the 
ERO number counts can be calculated, once their cluster- 
ing amplitude is known. 
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Fig. 12. The histogram shows the frequency distributions of 
the number of EROs (R - Ks>5 and Ks < 18.8) that have 
been recovered by sampling our survey with a SOFI field (25 
arcmin 2 ). The mean expected value is 10. The dotted curve 
represents the probability distribution expected in the poisso- 
nian case. 



To verify the consistency of this picture, we derived the 
distribution of the number of EROs (with R — Ks>5 and 
Ks <18.8, i.e. those in Fig. |J) that can be recovered in our 
area by sampling it with a field of view of 5'x5', which is 
the typical field of view of a near-infrared imager such as 
SOFI. In Fig. 12 the observed frequencies of the number 



of EROs recovered in this counts-in-cell analysis is plot- 
ted. As the mean expected number of EROs is about 10, 
the poissonian fluctuations would be a po isson ~3.2, while 
fluctuations with a =5.4 are actually observed. Applying 
Eq. (P, the measured clustering amplitude A — 0.013 im- 
plies a = 5.55, in excellent agreement with the measured 
a value. 

We also note that the distribution of the numbers of 
EROs in Fig. [l2] is not only asymmetric, but also very 
broad, ranging from A=0 to A=30. In 29% of the cases 
the number of EROs recovered is N <5, corresponding to 
a surface density half of the real one, while only in 19% of 
the cases the observed number is N >15. This shows that, 
on average, it is more probable to underestimate the real 
surface density of these objects. This is a clear property of 
the sky distribution that we observe, as the voids extend 
on a large fraction of the surveyed area. These results 
show how strong the effects of the field-to-field variations 
are in the estimate of the sky surface density of EROs. 
In this respect, it should be noted here that all previous 
estimates of the number density of high-z ellipticals were 



based on surveys made with small fields of view, typically 
ranging from 1 arcmin 2 in the case of the NICMOS HDF- 
S (Benitez et al. 1999) to 60 arcmin 2 in the case of Barger 
ct al. (2000). 

6.3. Implications for the evolution of elliptical galaxies 

The selection of galaxies with colors R — Ks > 5 can be 
used to search for elliptical galaxies at z > 0.9 (see Fig. 
|l3| ), and to study their evolution by comparing their ob- 
served surface densities with those expected from PLE or 
hierarchical models of massive galaxy evolution. In this 
respect, very discrepant results have been obtained so far, 
making the formation of spheroids one of the most con- 
troversial problems of galaxy evolution (see the Introduc- 
tion) . 

Our results on the ERO clustering clearly show that for 
such a comparison to be reliable, both a wide field survey 
(resulting in a large number of EROs) and a consistent 
estimate of their surface density fluctuations are neces- 
sary before reaching solid conclusions on the evolution of 
elliptical galaxies. 

In this section, with the main aim to show the effect 
of the increased uncertanties due to the clustering, a pre- 
liminary comparison is presented between the sky density 
of EROs observed in our survey (Table |) and the pre- 
dictions of an extreme PLE model similar to that used 
by Zepf (1997). In this model, ellipticals formed at Z{=5 
and their star formation rate (SFR) is characterized by 
an exponentially decaying burst with SFR oc exp(t/r), 
with r = 0.1 Gyr. Adopting the Markze et al. (1994) lo- 
cal luminosity function of ellipticals, and the Bruzual & 
Chariot (1997) models with solar metallicity and Salpeter 
IMF, the expected surface densities of passively evolving 
ellipticals with R — Ks > 6 was calculated for different 
limiting Ks magnitudes. 

Fig. [l4| shows the comparison between the expected 
and the observed densities of EROs with R — Ks > 6 (such 
a color threshold should select passively evolving galaxies 
at z 1.3). For each data point we show three different 
error bars, which are actually the region of confidence in 
the poissonian case (at la) and in the true (i.e. clustering 
corrected) case (at la and 2a). Such confidence regions 
have been estimated, following the prescriptions of Eq. (^|), 
by finding the range of values for the true average counts 
n for which the observed N would represent a deviation of 
the required number of a from the real density. In other 
words such ranges are defined from the two solutions of 
the equation: 



(n - A) 2 
n(l + nAC) 



(10) 



where a is the number of a considered. The amplitudes of 
the angular correlation function used for the R — Ks > 6 
EROs are those derived for the R — Ks > 5 EROs, which 
is likely to be a conservative assumption as the amplitudes 
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Fig. 13. The color-redshift relation for the extreme PLE 
model described in the test, computed for a fiat and an open 
universe (A = 0). The horizontal dashed lines are the thresh- 
olds adopted in this paper. The vertical line shows that z > 1 
ellipticals correspond to the R — Ks > 5.3 EROs. At z = 1 a 
Ks = 19 elliptical galaxy has L ~ 0.3-L* and L ~ 0.5L* in our 
flat and open models, respectively. 



Fig. 14. The observed integral number densities of EROs with 
R-Ks>6 (Table ||) are compared with the extreme PLE model 
described in the text. The three error bars shown here are, 
from left to right, the la poisson uncertainty and the 2a and 
la true uncertanties. The true uncertanties take into account 
the observed clustering of EROs. 



of redder samples should be higher, as suggested by Fig. 

^ in 

Fig. [L4| shows that the observed EROs densities are 
indeed lower than the predictions of this particular PLE 
model. However, even for the most deviant point, this PLE 
model can be rejected at only the 2.5a and 2.3<r level for 
f2o = 1 and = 0.1 respectively, if we use the "correct" 
error bars. Note that the different points plotted in that 
figure are not statistically independent because they are 
partially obtained with the same objects (they are cumu- 
lative values). 

It should be recalled here that the observed ERO den- 
sities plotted in Fig. 14 are an overestimate of the true 
density of the ellipticals because of the contamination by 
dusty starbursts (see Cimatti et al. 1998, 1999; Dey et al. 
1999; Smail et al. 1999) and by field low-mass stars. The 
fraction of dusty starbursts in complete ERO samples is 
not known yet, as discussed in the introduction, but our 
results show that they should not be the dominant pop- 
ulation. For instance, assuming that the fraction of dusty 
starbursts and low-mass stars is 20% and 10% of the ERO 
population respectively, this would decrease the observed 
densities plotted in Fig. [l4| accordingly, but it would also 
increase by a factor of 2 the clustering amplitudes of the 
high- z ellipticals (see Sect. |s~l| ), and hence the error bars 
related to those points. As a consequence, the statistical 



significance of the difference between data and model in 
this case would be only at the 2<r level. 

It is relevant to mention that the predictions of the 
PLE models depend very strongly on many parameters 
that have to be adopted a priori such as fioj the local 
LF of ellipticals (uncertain by up to a factor of 2), the 
redshift of formation Zf, the history of star formation, the 
metallicity, the IMF, the spectral synthesis models. For 
instance, even just a decrease of Zf, or a small residual 
star formation at z ~ 1.5 (Menanteau et al. 98, Jimenez 
et al. 99), would decrease the predicted numbers of EROs 
making them more consistent with our data. We therefore 
conclude that it seems premature to reject even extreme 
PLE models at a high level of statistical significance on 
the basis of these data. 

A preliminary comparison of our results can be made 
with some aspects of the hierarchical models of galaxy 
formation. First of all, our findings could qualitatively fit 
into the predictions of such models, where high-z ellipti- 
cals should be very clustered (Kauffmann et al. 1999) be- 
cause they are expected to be linked to the most massive 
dark matter haloes which are strongly clustered at high-z. 
The indication (marginally significant at ~ 2.7c level) of a 
decrease of the clustering amplitude of the EROs with the 



Ks magnitude (see Sect. 5.3), if mainly due to the mass 



of the galaxies, could also fit well in this framework be- 
cause smaller objects should be connected to smaller dark 
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matter haloes which are expected to be less correlated. 
On the other hand, our results seem to conflict with the 
predictions made by Kauffmann & Chariot (1998) on the 
fraction of -ftT-selected galaxies with K < 19. In fact, the 
fraction of galaxies observed to have color R — Ks > 5.3 
(which corresponds to the selection of z ^ 1 ellipticals) is 
about 7% of the total in our survey (see Table ||), to be 
compared with the 2-3% of z > 1 galaxies with K<19 
expected in the Kauffmann & Chariot (1998) hierarchi- 
cal model. This result on the fraction of z £ 1 galaxies in 
our Jf-selected sample broadly agree with the finding of 
Eisenhardt et al. (2000). 

7. Summary 

The main results of this work are: 

- We have presented a survey which covers 701 arcmin 2 
and is 85% complete to Ks < 18.8 over the whole area 
and to Ks < 19.2 over 447.5 arcmin 2 ; the R-band limit 
is R > 26.2 at the 3er level. 

- The observed galaxy counts are derived over the largest 
area so far published in the range of 18 < Ks < 19.2. 
Such counts are in excellent agreement with other pub- 
lished data. 

- The median R — Ks color of field galaxies increases by 
0.5 mags from Ks — 16.5 to Ks = 18, and it remains 
constant to Ks — 19.2. 

- A sample of 398 EROs has been selected. This sample 
is the largest published to date and is characterized by 
an area larger by about four times than previous sur- 
veys. The ERO counts and surface densities have been 
derived for several color thresholds and Ks limiting 
magnitudes. In particular, we find 0.67 ±0.03 (poisso- 
nian) EROs arcmin- 2 with R-Ks>5 and 0.10±0.01 
EROs aremkr 2 with R - K s > 6 at Ks < 19.2. 

- The surface density of EROs with R — Ks > 7 has 
been estimated for the first time to be of the order of 
~ 0.01 arcmin -2 at Ks ~ 19. 

- The angular correlation function of field galaxies, fit- 
ted with a fixed slope S — 0.8, has an amplitude 
A(l°) ~ 0.0015 at 18.5 < Ks < 19.2, in agreement 
with previous measurements. 

- For the first time, we detected the clustering of EROs, 
with an amplitude A(l°) ~ 0.015 for the objects with 
R - Ks > 5, in the range 18.5 < Ks < 19.2) which is 
about a factor of ten higher than that of field galaxies. 
The ERO two point correlations are very well fitted by 
a <5 = 0.8 power law. 

- The clustering amplitude of the galaxies increases with 
the R — Ks color threshold following the relation 
log A oc 0.43(i? - Ks), for 3 < R - Ks < 5.7 at 
Ks < 18.8. 

- The strong clustering of EROs is shown to be a direct 
evidence that a large fraction of these objects are in- 
deed high-^ ellipticals. Our result is therefore the first 



detection of the large scale structure traced by the el- 
liptical galaxies at z ~ 1. 

- The ERO clustering explains the conflicting results ob- 
tained so far on the density of high- z ellipticals in terms 
of strong field-to-field variations affecting the surveys 
based on small fields of view (e.g. 5x5 arcmin). 

- Taking into account the clustering of EROs, even the 
predictions of extreme PLE models for the comoving 
density of high- z ellipticals cannot be rejected at much 
more than 2a significance level. 
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